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Synthetic quantum orbital materials, such as 
the recent double-well optical lattices loaded with 
s and p orbital atoms jl|, 0], open an avenue to 
exploit symmetries beyond natural crystals. Ex- 
otic superfluid states were reported for bosons. 
Here, we unveil a topological phase of interact- 
ing fermions on a two-leg ladder of unequal par- 
ity orbitals, derived from the experimentally re- 
alized double-well lattices by dimension reduc- 
tion. Z2 topological invariant originates simply 
from the staggered phases of sp-orbital quantum 
tunneling, requiring none of the previously known 
mechanisms such as spin-orbit coupling or artifi- 
cial gauge field. Another unique feature is that 
upon crossing over to two dimensions with cou- 
pled ladders, the edge modes from each ladder 
form a parity-protected flat band at zero energy, 
opening the route to strongly correlated states 
controlled by interactions. Experimental signa- 
tures are found in density correlations and phase 
transitions to trivial band and Mott insulators. 

The uneven double- well potential formed by laser light 
has been developed into a powerful tool for quantum 
gases by numerous groups [Irl6| • For instance, controls of 
atoms on the s- and p-orbitals of the checkerboard [l[ and 
hexagonal @ optical lattices have been demonstrated, 
and correlation between these orbitals tends to give ex- 
otic quantum states @, 0). Motivated by these devel- 
opments, we consider a lattice of uneven double- wells 
where fermionic atoms are loaded up to the s- and p- 
orbital levels of the shallow and deep wells respectively, 
as shown in Fig. [1] The spatial symmetry of the orbital 
wave function dictates the complex hopping amplitudes 
between nearby sites. Under certain circumstances, as 
for the uneven double wells, the orbital hopping pattern 
is sufficient for producing topologically nontrivial band 
structures @. 

We will first focus on a one-dimensional ladder sys- 
tem illustrated in Fig. [Tp and \Tjp. This corresponds to 
the quasi-one dimensional limit of a standard double- well 
optical lattice, with the optical potential given by 



V(x,y) = V x sin 2 (kx) + Vism 2 (ky) + V 2 sm 2 (2ky +-^). 
This optical lattice has a double well structure in the 



laser 




X A 




-TV 



7T 



*email: w.vinccnt.liu@gmail.com 



FIG. 1: The uneven sp-orbital ladder system made from a 
two-dimensional double-well optical lattice through dimen- 
sion reduction, a, An optical lattice of uneven sub- wells (light 
and dark blue), with parameters V x /Vi = 0.3, V2/V1 = 1 and 
<j> = 0.67T, develops high barriers (red ridges) in the y direc- 
tion, slicing the lattice into dynamically decoupled uneven 
two-leg ladders, b and c, Schematic side and top views, re- 
spectively, of the ladder illustrate tunneling (t's) of fermions 
prepared in the degenerate s and p x levels, d, Topological 
winding of Hamiltonian across the Brillouin zone. 



y-direction. For V1.2 S> V x , there is a large tunneling 
barrier between double wells in the y-direction, so in low 
energy physics the two dimensional system decouples into 
an array of dynamically isolated two-leg ladders of A and 
B sub- wells (Fig. []}, with each ladder extending in the 
x-direction. The relative well depth of the two legs is con- 
trolled by the phase 4> and further by the ratio V2/V1. We 
will focus on a situation, similar to the setup in the ex- 
periment [l[ , where the s-orbital of leg A has roughly the 
same energy as the p^-orbital of leg B (other p-orbitals 
have much higher energy). For example, one can choose 
Vi = 40E R , V 2 = 20E R , V x = 3.8E R and = 0.9tt 
in experiments, where E R is the recoil energy h 2 k 2 /2m, 
with m the mass of atom and k the wave number of the 
laser. Such a setup will give the A (B) wells a depth 
2.7E R (8.1Er). The lattice constant a = n/k will be 
set as the length unit in this paper. We now consider a 
single species of fermions occupying these orbitals, with 
the low-lying s-orbital of leg B completely filled (altcrna- 
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tively fcrmions can be directly loaded into the p^-orbital 
of leg B by techniques developed for bosons in Ref. [H,|I|]). 
The Hamiltonian of the sp-orbital ladder is then given 

by 



by introducing auxiliary operators, <j>±(j) 



c 



Y.j l>C:Cr (1) 



where C 



the 



fermion creation operators for the s- and p^-orbitals on 
the A and B leg respectively. The relative sign of the 
hopping amplitudes is fixed by parity symmetry of the 
s and p x orbital wave functions. As depicted in Fig. \Vp, 
the hopping pattern plays a central role in producing a 
Z 2 topological phase. With a proper global gauge choice, 
t s , t p and t sp are all positive. The rung index j runs from 
to L— 1 with L the system size. We consider half filling 
(one particle per unit cell) , for which the chemical poten- 
tial /I = 0, and the Hamiltonian is particlc-hole symmet- 
ric under transformation Cj — > (— l) J Cj. This sp-orbital 
ladder contains a staggered hopping, which is analogous 
to the Su-Schricffer-Heeger model for the dimerized poly- 
acetylene, an electron-conducting polymer [9|. Its physics 
is also connected to the more familiar frustrated ladder 
with magnetic flux [To| , but the sp-orbital ladder appears 
much easier to realize experimentally. 

In the momentum space the Hamiltonian takes a sim- 
ple and suggestive form, 



H(k) = ho(k)I + h(k)-a, 



(2) 

= 2t sp sin(fc) 
is the unit 

matrix, a x , a v and a z are Pauli matrices in the two- 
dimensional orbital space. The energy spectrum consists 
of two branches, 



where ho(k) = (t p —t s ) cos(fc), h x ■ 
and h z (k) = — (t p + i s )cos(fe). 



0, h y (k) 
Here, I 



E ± (k) = ho{k) ± yjhl(k)+hl(k), 



(3) 



with a band gap E g = min(2i p + 2t s ,4t sp ), which closes 
at either t sp = or t s + t p = 0. An interesting limit that 
highlights the nontrivial band structure of our model is 
that when t p = t s = t sp , the two bands are both com- 
pletely flat. To visualize the topological properties of 
the band structure, one notices that as k is varied from 
— 7r through to +7r, crossing the entire Brillouin zone, 
the direction of the h(k) vector winds an angle of 2ir 
(Fig. [T}1). The corresponding Berry phase is half of the 
angle, 7 = it. For a one-dimensional Hamiltonian with 
particle-hole symmetry, the Berry phase is quantized fll| , 
7 mod 27r = or n, which defines the trivial and Z2 
topological insulator, respectively. The sp-orbital ladder 
exhibits 7 = 7r for t sp 7^ 0, so it is a Z2 topological insu- 
lator at half filling. This is verified by explicit calculation 
of the Berry phase (see Supplementary Information). 

The nontrivial topology of the ladder system also man- 
ifests in existence of edge states. It is easiest to show the 
edge states in the flat band limit, t s = t p = t sp = t, 



a s (j)]/y/2. Then the Hamiltonian only contains coupling 
between </> + and 0_ of nearest neighbors, but not among 
the </)+ (or <j)J) modes themselves, 



H ^2tJ2 i 4>l(j)4> + (j + l)+h.c. 



(4) 



Immediately, one sees that the operators </>+(0) and 
4>-{L — 1) at the left and right end are each dynami- 
cally isolated from the bulk, and do not couple to the 
rest of the system (Fig. Wt*)- These loners describe the 
two edge states at zero energy. 

It is instructive to map the sp-orbital ladder to two 
decoupled Majorana chains [l2|. For each rung j, we 
introduce 4 Majorana fermion operators, if)\± = (4>± + 
h.c.) and ip2± = (*</>=f + h.c). Then the Hamiltonian is 
transformed into 

Ho -> itJ2j h/>2+0>i+(j + 1) - Vi-(j)V>2-(j + 1)] , 

(5) 

for which four unpaired Majorana fermions ?/>i+(0), 
ip2-(0), tp2+(L — l) and tpi-(L— 1) emerge. The four un- 
paired Majorana fermions manifest a Z2 (Z4) degenerate 
ground state in the canonical (grand-canonical) ensemble 
at half filling. Note that the Majorana number as defined 
in Ref. [T3] is 1 for our double Majorana chains. While 
on general grounds coupling between the four unpaired 
Majorana fermions is allowed, local coupling, e.g., in the 

form of ieVi+ (0)^2- (0) = — 2e ^Cq [<tq + o x \ Co^j + const, 
is however prohibited in our system due to the required 
particle-hole symmetry. 

Once the particle-hole symmetry is broken, topologi- 
cally protected Majorana fermions with Majorana num- 
ber — 1 [12] can be realized on the s p- ladder using schemes 
similar to those proposed in Ref. |13l - fl5| . e.g., by induc- 
ing weak pairing of the form ■ Aa s (j)a p (j) + h.c. In 
the parameter regime 2t s < \u\ < 2t p , the Majorana 
number is —1 and the resulting Majorana zero modes 
are topologically protected. The topologically protected 
Majorana state is a pro mising candidate for topological 
quantum computing [l6l . [13] ■ Compared to previous pro- 
posals p^Hla ]. the present ladder system does not re- 
quire spin-orbit coupling. 

Away from the flat band limit, in general the wave 
functions of the edge states are no longer confined strictly 
at j = or L — 1, but instead decay exponentially 
into the bulk with a characteristic length scale £ = 
2 / 1 °g(l(v / ** i p + t sp)/(\/t s t p - t sp )\). The analytical ex- 
pression for the edge state wave functions in the general 
case are discussed in Supplementary Information. At the 



critical point t s 



0, the bulk gap closes and £ — > 00. 



The existence of zero energy edge states is also confirmed 
by numerical calculation as shown in Fig. [2}d and [2):. 

For a finite ladder of length L with open boundary 
condition and populated by L fermions (half filling), 
L — 1 fermions will occupy the valence band (bulk states) 
and one fermion will occupy the edge states (Fig. [2]). 
Since the two edge states are degenerate, the ground 
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FIG. 2: Bulk and edge eigenstates of the sp-orbital ladder, 
a. A pictorial representation of the simplified Hamiltonian 



in the flat band limit t s = 



t ap showing the emergence 



gap scales as exp(— aL) as L — > oo fl8j . 

An interesting topological phase transition to a triv- 
ial insulator can be tuned to occur when rotating the 
atoms on the individual sites, for ex amp le, by applying 
the technique demonstrated in Ref. |19j. Such an indi- 
vidual site rotation amounts to addition of an imaginary 
transverse (along y) tunneling between s- and p^-orbitals 
in our Hamiltonian, SH = A w ■ Cja y Cj. This term 
preserves particle-hole symmetry but breaks both parity 
and time-reversal symmetries. The total Hamiltonian in 
the momentum space now reads "H'(fc) = T-L(k) + A y a y . 
The quantized Berry phase remains at ir for A y smaller 
than the critical value A y = 2t sp , so the topological 
insulator phase is stable against finite A y . This fur- 
ther shows that the Z 2 topological phase is protected 
by particle-hole symmetry. For A y greater than A y , 
Berry phase vanishes and the system becomes a triv- 
ial band insulator. At the critical point the band gap 
closes. Apart from the Berry phase, the topological dis- 
tinction between 7i'(k) and %{k) can also be seen from 
a gapped interpolation (Tl| as discussed in Supplemen- 
tary Information. Besides probing the half charges on 
the boundaries, another signature for the critical point 
of the topological phase t ransition i s the local density 
fluctuation, dp = j E i V (pti)p(j))- $P 1S l/v^ when 
A y = 0, independent of other parameters t s , t p and t sp , 



of isolated edge modes. The definition of the <fi± operator 
is given in the main text. b. The eigen energy of a ladder 
with finite length L — 12 showing two degenerate zero energy 
states inside the gap. c. The probability distribution of the 
in-gap states (Eq. ©) for varying strengths of inter-orbital 
interaction U sp . The in-gap states are shown localized on the and decreases monotonically with increasing A y . The 
edges and survive against finite interaction. In b and c, we peaks of dSp 2 /dA y reveal the critical points (Fig. [3)d) 
choose t s — t p — 2t sp (taken as the energy unit). and provide a reliable tool of detecting the topological 

phase transition in experiments. 

It is feasible to prepare the ladder with phase separa- 
tion: e.g., a topological insulator on the left half but a 
trivial insulator on the right half. This can be achieved 
by rotating the lattice sites on half of the ladder only. 
The system is now described by 



state has a double (Z 2 ) degeneracy. The edge state is 
a fractionalized object carrying half charge (cold atoms 
are charge neutral, here charge refers to the number of 
atoms). This becomes apparent if we break the particle- 
hole symmetry by going infinitcsimally away from half 
filling, e.g., tuning chemical potential p = + . Then, the 
valence bands and the two edge states will be occupied. 
With a charge density distribution on top of the half 
filled background defined as p(j) = C^Cj — 1, one finds 

Eo<j<d(/>(i))U=o+ = Ei-d-x<j'<L(^C/'))U=o+ = h 
where d satisfies ( « d « L ( c -g-, take d = 5£). A char- 
acteristic feature of the Z2 topological insulator (with the 
number of atoms fixed) is the topological anti-correlation 
of the charge at the boundaries, 



H v = H + ^'£ j [l- cosr?(j)] C}a v C. 



(6) 



lim J2(p(j)p{j')) 



In the flat band limit, £ — > 0, the edge states are well 
localized at the two ends of the ladder. The topological 
anti-correlation simplifies as (p(0)p(L — 1)) = —i, and 
the half charge is also well localized, i.e., (p(0)) M=0 + = 
(p(L — 1)) M=0 + = \- Since the edge states are well iso- 
lated from the bulk states by an energy gap, they are 
stable against local Gaussian fluctuations. The coupling 
between the two edge states vanishes in the thermody- 
namic limit (L — > 00), because the hybridization induced 



with a field configuration rj(j) which satisfies the bound- 
ary conditions rj(j = —00) = and rj{j = +00) = 7T, 
and A y > A y . The charge distribution induced by the 
domain wall (the phase boundary) is calculated both 
numerically and from effective field theory shown in 
the Supplementary Information. Both approaches cross- 
verify that the domain wall carries half charge (Fig. [3^) . 
The half charge can be detected by the single site imaging 
technique in experiments [20l . l2l| . 

We further examine the stability of the topological 
phase and its quantum phase transitions in the presence 
of interaction using exact diagonalization. For single- 
species fermions on the sp-orbital ladder, the leading in- 
teraction term is the on-site repulsion between different 
orbitals, 



H- m t — U sp 



n sU) - ^ 



n p(j) - \ 



(7) 



We compute the fidelity metric g as function of the inter- 
action strength U sp (see the Methods section). A peak 
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FIG. 3: Phase transition between the topological and the triv- 
ial band insulator, a. A domain wall between a topological 
insulator (t s = t p = t sp , A y — 0, left) and a trivial insulator 
(t s = t p = t ap = 0, right). The circle represents the delo- 
calized fermion shared by two neighboring rungs as depicted 
in Fig. Oi, whereas the ellipse represents localized fermion 
without hopping. The additional charge i in the middle is 
the fractional charge carried on the domain wall. b. The 
derivative of density fluctuation, — It develops sharp 

peaks, measurable in experiments, along the line of topologi- 
cal critical points. 
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FIG. 4: Transition from a topological insulator (TI) phase 
into a Mott insulator with ferro-orbital (FO) order with in- 
creasing interaction. Top panel shows the fidelity metric g 
and the ferro-orbital order parameter X 3p . Bottom panel 
shows the particle/hole chemical potential (n P /nh). The 
finite charge gap — /Xh in the bulk calculated with peri- 
odic boundary condition (dashed lines) comparing with the 
vanishing gap with open boundary condition (solid lines) in- 
dicates in-gap states on the edge. The length L is 12, and 
t s = tp = 2t ap (taken as the energy unit) in this plot. 



in the fidelity metric indicates a quantum phase transi- 
tion [22j ■ Our numerical results (Fig. [3J show that the 
topological phase is stable for U sp < U^ p , with robust 
in-gap (zero-energy) states [22| localized on the edges 
(Fig. For stronger interaction the ladder undergoes 
a quantum phase transition to a Mott insulator phase, 
which exhibits ferro-orbital order with order parameter 
defined as ^X sp = {C'-cx x Cj). 

Finally we point out a remarkable difference from the 
Su-Schricffer-Heeger model. The zero-energy states in 
the present sp-ladder survive even when the system is 
extended to two dimensions with finite inter-ladder cou- 
pling (e.g., by reducing V\ t 2 relative to V x in the setup of 
Fig. Qli). The zero modes of individual ladder morph into 
a flat band with double degeneracy (Fig. [5]). The lack of 
dispersion in the y direction is related to the inter-ladder 
hopping pattern, which does not directly couple the edge 
states but only s- and p-orbitals on different rungs. The 
unexpected flat band in 2D is an exact consequence of 
the 7>orbital parity and hence is protected by symmetry. 
The flat dispersion can be rigorously proved using an uni- 
tary transformation and arguments based on the quanti- 
zation of Berry phase in the Supplementary Information. 
The flat band is reminiscent of that at the zigzag edge 
of graphene, but with the difference that the present flat 
band is protected by the parity of the orbital wavcfunc- 
tions. The diverging density of states associated with the 
flat band provides a fertile ground for interaction-driven 




FIG. 5: The energy spectra of the two dimensional system 
of coupled ladders with t s = t p = 2t 3p (taken as the energy 
unit here), and length L = 200. An open (periodic) boundary 
condition is applied in the x (y) direction, a and b show 
the spectra with the small and large inter-ladder coupling, 
t'sp = tsp/5 and t' sp = t sp , respectively. A flat band (red line) 
at zero energy with double degeneracy generically appear for 
^ £ s p <C t sp - 



many body instabilities. Future work will tell whether 
strongly correlated topological states exist in such two- 
dimensional interacting system. 
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Methods 

Imaginary transverse tunneling. By rotating in- 
dividual lattice sites the induced bare coupling term 
is ^2jQL z (j) where the angular momentum operator 

L z = —i(a* p a Px — a Px a Py ). This term couples the p x 
to py-orbitals of the B-leg. One can tune the rotating 
frequency to match f2 with the transverse tunneling t sy 
from the pj,-orbital of the B-leg to the nearby s-orbital on 
the A-\eg. Despite the large energy band gap (e y ) which 
separates the p y -orbitals from the degenerate s and p x 
orbitals (bare in mind the s and p x orbitals are from dif- 
ferent legs of the ladder) , the low energy effective Hamil- 
tonian receives a standard 2nd-order effect from virtual 
processes, in which a particle jumps from a p^-orbital to 
the on-site Py-orbital and then to the nearby s-orbital. 
The correction is given by —C^a y Cj, which makes an 
imaginary transverse tunneling between nearby s- and 
p^-orbitals. 

Interaction effects. To characterize the stability 
of the topological phase against the inter-orbital in- 
teraction Hi n t (see Eq. (JT))), we use the exact di- 
agonalization method to calculate the fidelity met- 
ric g ee 2[l~\^l(U sp )\^° L (U sp + 5U sp ))\]/L(8U sp ) 2 : 
where \iP°l{U sp )) is the ground state wave function of 
the Hamiltonian H = Hq + H- m t for a finite chain of 
length L. In presence of interaction, the edge states sur- 
vive as in- gap states (zero energy single particle/hole 
excitations) [22(. The energy of single particle (hole) 
excitation is defined as fi p (fj,h) by a p = El+i — El 
(fih = El — El~i), where En is the ground state en- 
ergy of the ladder loaded with N fermions. The spatial 
distribution of the in-gap states is defined as the density 
profile (Arij) of a hole created out of the ground state, 
which is 

An, = <V£|p(jM> - Wl-MjM-i), (8) 



where \tp N ) is the ground state with TV fermions. The 
density profile Arij is found to be localized on the edges 
when U sp <C Ug p , and to delocalize when approaching the 
critical point and finally disappear (Fig. [5fc) . The Mott 
state appearing at the strong coupling regime has a ferro- 
orbital order (6 sp (j)), withO sp (j) = C^azCj (Fig.EJl. In 
our numerical calculation of finite system size, the cor- 
relation matrix [C] Ji:)2 = (O sp (ji)O sp (j2)} is calculated 
and the strength of the ferro-orbital order A sp is defined 
as the maximum eigenvalue of [C]/L, extrapolated to the 
thermodynamic limit. 
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Supplementary Information 

S-l. TOPOLOGICAL INDEX 



The topological nature of the sp-orbital ladder can be understood in terms of the winding number of the Hamiltonian 
in the momentum space. Given a Hamiltonian 



H(k) = h (k)I + h(k)-a, 
with h x = 0, the winding number is defined as 

/dk 
—e vv -h~ 1 (k)d k h v -, 



(SI) 



(S2) 



-zy 



1. This winding number is 1 for the sp-orbital ladder in the topological insulator 



where h = -M- and e vz 
phase (Fig. [Til). 

The Hamiltonian in equation (|S1|) can be written as T~L(k) = U(l + a z )U\ where U — exp (icr x 8(k)/2), with 8(k) 


1 



defined by 



" cos(0) " 




' k(k) ' 


sin(0) 




. hy( k ) - 



The eigenvector of the lower branch is given as u(k) = e t9 ( fc )/ 2 £/ 



The Berry 
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phase is given by 



7 = <j)dkiu^dkU 



dk 

= - <£ dkd k 
= -Wit. 



-ia !c e/2 e i6/2Q k ^ e ia x 0/2 e -i9/2 



dk- 



dk sin(0) d k cos(0) 



cos(6») 



sin(0) 



(S3) 

With the eigenvector u(k) multiplied by a phase factor e l ^ k ^ in which § dkdk<t>{k) = 2mr, the Berry phase changes 
2mr. This means the Berry phase 7 = Wn mod 2tt. The winding number W being even and odd defines two classes 
of topological states in one dimension. 



S-2. THE WAVE FUNCTIONS OF EDGE STATES 



To see the existence and the property of of the edge states away from the flat band limit it is useful to rewrite H 
in the following form 



Ho = Ej *} 



with 



*-(?*) 



t a 



<I> 



h.c 



C h 



(S4) 



(S5) 



and t\ = (tp — t s ), ti = yjt s t p — t sp and £3 = yjt s t p +t sp . The wave functions wi/ rM (j) of the edge states are introduced 
by h /r = j2ju= s ,p w i/r^(j)av(.j) = 12j a =± ^i/r,a(j)^a(j)i where b l/r is the fermion operator of the left/right edge 
state. The wave function of the left edge state bi is given as 



CXp ; 

0, 



if j is even; 
otherwise, 



t2+t 3 



exp(-j/0 , if j is even; 



(S6) 



with the decay width £ 



0, otherwise, 
The wave function of the right edge state can be constructed by performing 



log(|t 3 /t 2 |) 

parity transformation of the left edge state. As long as t 2 vanishes, the width £ — >• and the edge states are completely 
confined on the ends of the ladder. With finite t 2 the edge states deconfine and decays exponentially. The critical 
point is £2/^3 = 1 (t S p = 0), for which £ — > 00. At this point the bulk gap closes (Eq. Q). The analysis presented 
here is confirmed by numerical calculations. 



S-3. MAPPING TO A MODEL OF TWO MAJORANA FERMI CHAINS AND DECONFINED 

EXCITATIONS 



In the flat band limit t s — t p — t sp = t (or equivalcntly t\ = t 2 = 0) we can introduce particle-hole mixed 
operators 4>^/\,{j) = [a p (j)±a\(j)]/ \/2~. The sp-orbital ladder maps to two decoupled Kitaev's p-wave superconducting 
chains lla. 



H -> tJ2 {4(7')^(j + 1) + P^UMHJ + 1) + h.c.} 



(S7) 



with P-i 



It is instructive to rewrite the fermion operators <f>± defined in Sec. IS-21 in terms of Majorana fermion operators as 
<j>± = 5 (ipi± — ^2^), with ipip — ip\ v , and {V^pi V^'p'} = 28ppi8u'- The resulting Hamiltonian in Eq. (|S4[) reads 



^0 = ^[^l + (j),V2-(j)^l-(j)>2+(.?)] 



- 


-h 





-*2" 


*1 





£2 








-*3 




















^i+(i + i) 
^2-(i + i) 
^i-(i + i) 
^2+(i + i) 



(S8) 
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In the flat band limit, the Hamiltonian maps to two decoupled Majorana Fermi chains [I2j . 

#o -> *f Ej + 1) - 4>i-U)ih-(j + 1)] ■ (S9) 

The unpaired Majorana fermion operators are - 0i+(O), i/) 2 _(0), ip2+(L — 1) and %pi-(L — 1). These Majorana fermion 
operators give a many-body ground state manifold of Z4 degeneracy for fi = in the grand canonical ensemble fl2l ]. 
Since [<f>\. (O)0+(O), if] = 0, [</>!(£ - - l),ifj = and (O)0+(O), 01(1, - 1)<^_(X - 1)] = 0, one can label 

the degenerate ground states by operators 0+(O)</>+(O) and <f>_(L — l)<fi_(L — 1) in such a way as follows, 

|G ): ^(0)^ + (0)|G )=Oand^(i-l)^-(i-l)|Go)=0, (S10) 

\G+) = 4(0)|G >, (Sll) 

|G_) = ^_(L-l)|Go>, (S12) 

|Gf+_) = ^(0)^L(L-l)|G o ). (S13) 

(The states \G+) and |G_) have the same number of fermions, providing the Z 2 degeneracy at half filling in the 
canonical ensemble.) In the Z4 degenerate ground state manifold the operators <^+(0) and </>_(L — 1) read as 

cf> + (0) = \G o ){G + \ + \G.){G + .\ (S14) 
ML-1) = |Go)(G_| - |G+)(G+_|. (S15) 

To demonstrate the deconfined Dirac fermion excitations explicitly in this ground state manifold, we can label the 
states by d\d\ and d\d 2 , where d\ and d 2 are fractionalized Dirac fermion operators d\ = (-0i+(O) — iip 2+ {L — l))/2 
and d 2 = (ip2-(0) ~ iipi-(L — l))/2, respectively. The relabeled states are defined by 

|G >: dld^Ga) = and d\d 2 \G Q ) = 0, (S16) 

|G 1 )=4|G >, (S17) 

|G 2 )=4|G ), (S18) 

|G 12 ) =44|G ). (S19) 

These states embedding fractionalized Dirac fermions are actually given by \Gq) = (|G+) — |G_))/\/2, |Gi) = 
(|G ) - |G+_))/V2, |G 2 ) = -(|G ) + |G+_»/\/2 and |G 12 ) = -(|G+) + \G-))/V2. A transition from |G ) to |Gr) 
creates one fractionalized Dirac fermion excitation. 

Away from the flat band limit the unpaired Majorana fermion operators tp are defined by [ip, Ho] = 0. The four 
unpaired Majorana fermion operators are ip\ = bi + bt, ip2 = i(pi — b\), tp^ = b r + b\ and ip4 = i{f> T — b\), where the 
Dirac fermion operators bi and b r are defined in Sec. IS-21 



S-4. THE GAPPED INTERPOLATION. 



An interpolation between %{k) and 'H'(k) is defined as h(k,<p) = W(k) + -f- [1 — cos (ip)] a y + A x sm(tp)cr x = 

do(k) + d(k, tp) ■ f?, with h(k, 0) = Ji(k) and h(k, n) = H'(k). The eigenvalues of this interpolation exhibit a finite gap 

for any (k,(p), given a sufficiently large A^. The Chern number C\ defined by C\ = J dkdipd ■ x is f for 

A y > Ay, and is for < A y < A y . Since the Chern parity ( — l) Cl is well defined for gapped interpolations with 
fixing h(k,0) and h(k,ir), C\ being odd tells H(k) is topologically distinguishable from H'(k) [ill ]. 



S-5. EFFECTIVE FIELD THEORY OF THE DOMAIN WALL 



To calculate the charge modulation induced by a domain wall connecting a topological insulator and a trivial 
insulator, we couple the ladder system to an auxiliary charge U(l) gauge field (A t ,A x ) and construct a low energy 
effective field theory, which is 1 + 1 dimensional reduced Chern-Simons Field theory (llj . A r\ dependent model 
Hamiltonian describing the domain wall is introduced 

H v = ^o + ^^Il-cos^Gj^Gj + A^sin^Gt^Gj (S20) 
3 

= ^ C k H v (k)C k . (S21) 

k 
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Here the a x term is introduced to guarantee a finite energy gap (A e ) of H v . We split the r\ field into two parts as 
V = Vo + ^v{.h T )- The field t] satisfies the boundary condition — 7y (— oo,r) = 0, 770(00, r) = ir. For example it can 

be set as 770 (j, r) = arctan + f j where ld w is the characteristic width of the domain wall. To proceed we assume 

Idw is much larger than the microscopic length scale (~ l/ z ) of the fermion system. The field Srj(J,T) satisfies a 

periodic boundary condition. In the following derivation we will treat 770 as a quasi-static field since it varies slowly 
in space. (One can think that we are deriving a local effective field theory for a subsystem in which the variation of 
770 is negligible.) 

With the model Hamiltonian H v coupled to the gauge field, its Lagrangian is given by 

L[C\C]= / dr£(r) 
Jo 

C = E [c]D T Ci + C}T JJ+1 C J+1 + C} +l T ]+1] C 3 + ^ [1 - cos(77)] f>,r ; + A, sinfaJC^C,-} 



(S22) 



with 



D T = d T + iA T (j,r), 



3 3 + 1 



^sp 

-t 



It is readily verified that the Lagrangian has a t/(l) gauge symmetry defined by 

A T -> A T + 9 T ??(i, r) 

A x (j + 1/2, r) -> A^j + 1/2, r) + + 1, r) - r)] 
C7(j,-r)->CC7,r)e-^». 



(S23) 



(S24) 



In the continuum limit, A x transforms as A x (x, r) — > A x (a;, r) + d x ip(x, r). In the following we shall write fc as fc x , 770 
as k y and ^77 as A y just to make the equations more compact. 

In the momentum-frequency (K = {k x ,iuj)) space the Lagrangian reads 



KK 1 



®KK' = (- iuj + Uri(k x ))6 K K> 



5V K k> = ^^{K - K')T V ( 



k x + k' x 



(S25) 



where F T = il, T a=x ^ y = ^? ■ For convenience, we further introduce G(K) = [—iuj + H v (k)] 1 
The effective action of the gauge field A v is defined as 

S cS [A, y ] = - In J D[C\C] cxp {-L[C\C- A v ]) , 

which can be calculated order by order as follows. 

S e s[A v ] = - hidet[X>] = -Tr InD = -Trln [p (0) + SV 



-TrlnX> (0) - Tr 



1. 



-Tr 



(S26) 



(S27) 



Since we are calculating the charge/current induced by the domain wall configuration — 770 field, we only keep the 
quadratic terms and thus the effective action is given by 

2 



= \ E MQ)M-Q)^J2 G ( K + Q/2)T^K)G{K - Q/2)T V {K). 



(S28) 



A" 
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Although the theory is complicated, it can be simplified because we are only interested in a low energy theory. The 
generic form of S e g at low energy is 

Seff = I3LJ2MQ)M-Q) + iqoW^) , (S29) 

K 

with Q = (q,iqo)- The massive terms are prohibited due to gauge invariance. By comparing the above equations we 
get 

= ^Q^o^ Tr ^ G(/ ' +0/2)r ^ (X)G( ^ _Q/2)r,y(X) ' (S30) 
KMV = ^^Tq TT ^ GiK + Q/2)T ^ K)G{K ^ Q/2)TAK) - (S31) 



K 



Then we have 



W" = -^— r TvY,[d l ^G{K)]T^K)G{K)T y {K)-^^ v , (S32) 

KT = ■^ I TtY J [dkG(K)\T IJl {K)G{K)T v {K)- l x^v. (S33) 

It is clear that W and JC are anti-symmetric. From Ref. [TlJ] , W xy = § § 4^ = i7^7( ? /) where O feai?) (7(77)) 
is the Berry curvature (Berry phase) of the Hamiltonian H n . Other terms of W — W Tl/ — vanish because G(K) and 
diujG(K) commute. The commuting relation can be verified by choosing the eigen-basis of H n (k). Now we have fully 
established the frequency part of the effective action, which is 

PL A p (K)A V (K)iuW u = 2(3LY / W xy A x (-K)(-iu)A y (K) (S34) 

K K 

In the real space (continuum limit) this term reads 2 J dxdrW yx {A y d T A x ). From gauge invariance, it is readily proved 
that —iJC Ty = W xy . Thus the effective action is given by 

S eB [A^] = 2W yx J dxdT(A y d T A x - A y d x A T ). (S35) 

The real time action reads 

S eS [A„] = -2iW yx J dxdt(A y d t A x - A y d x A t ) 

= f dxdt(A x dtJ] — A t d x i r i)-^—d ri "f(ri). (S36) 

We have replaced d v 8rj by d v r\ because d v r] w d v 8r\. This action can actually be written down directly from the famous 
2D Chern-Simons field theory by a dimension reduction procedure (Tlj . The linear response of charge/current is given 
b y 3" = Ta^- The cnar S e carried by the domain wall is given by Q = / fjg = J dxd^d^ir)) = - f ^^7(77). 
The charge Q is 1/2 with A x > 0. The microscopic details of the domain wall can change Q by 1, so we conclude the 
charge carried by the domain wall is Q = 5 mod 1. 



S-6. LOCAL DENSITY FLUCTUATION OF TOPOLOGICAL (NON)-TRIVIAL INSULATORS 

In this section we derive an analytic formula for the density fluctuation, which can be measured in cold atom 

experiments. The density fluctuation in cold atom experiments measures Sp = y~] j (^(cjCj — 1^ 

First two flat-band limits — (1) (A y = 0, t s = t p = t sp 7^ 0) and (2) (A y 0, t s = t p — t sp = 0) arc explored. The 
two limits give different Berry phases (7 = ir for Limit (1) and 7 = for Limit (2)). Limits (1) and (2) give non-trivial 
and trivial insulators respectively. For limit (1) fermions live on the bonds (one fermion per bond). The ground state 
can be written as 

l G > = IT 7| - Mi + 1)) 10), (S37) 
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with |0) the vacuum state with no particles. ((CjCj) 2 ) is given as follows 



= E {^'(<^0')<Mj)> - s^-a'i&MlWMfaU))} 

cx,<x' 



1 + 2 x 



(S38) 



Then we have 5p = For Limit (2) there is no hopping between different unit cells and fermions localize on each 

unit cell. Fermions cannot tunnel between different unit cells, so the local density fluctuation vanishes, i.e., bp = 0. 

For the generic case, the periodic boundary condition is adopted and the calculation is performed in the momentum 

space. In momentum space, the Fourier transformed operators are defined as a (k) = 4^J2j a j e ~ ikj - The density 

fluctuation, which is a bulk property, does not depend on the boundary condition in the thermodynamic limit. The 

r- -\T 

fermion operators of the eigen-modes (labeled by k) are introduced by b\(k),bi(k) = WC k - t/4- nere means the 
upper/lower band. The unitary matrix U is defined in Sec. IS- 1 1 The ground state of the fermionic ladder at half 



filling is \G) = Uk bl{k)\0)- The calculation of {{C]C 3 ) 2 ) is as follows 



L 2 ' 



\i E E ^^ +fe *-^- fci ^(4(fc 1 )5 t/ (fc 2 )a|:,(fc 3 )5^(fc 4 )) 

^EE{^(fc)a4fc)>(4'(^^(fcO) + (4(fc)^(fc))(a,(fc')4'(fc'))} 

kk' 

^E I (bl(k)h(k)b\(k%(k')) + e ^ 1 ( fc )^( fc >L(fc')^(fc)>[f/(fc)]L 1 ^(fc')]- 2 ^(fc')]:'.s 3 [^(fc)],' S4 

kk' I VV 1 ,S\S2S3S4 



It 



(S39) 



The term -k ~£ kk ' \[ u \ k )U(k')]^\ 2 simplifies as 



± 2 Y\PHk)u(k% t \ 2 = ^Y\h 



-ia x 9(k)/2 icr x e(fe')/2 



it 



1 v-- 2 f9(k) - 9(k') 
^E snr 



dk 
2^ 



cos{9(k)) - 



dk 
2^ 



sin(0(fc)) 



(S40) 



Due to the particle-hole symmetry h z (k) = —h z {-K — k) and that h(k) = h(ir — k), and we thus have cos($(fc)) = 
— cos(6>(7r — fc)), so that § dk cos(9(k)) = 0. Then the density fluctuation of a particle-hole symmetric insulator is 
given by 



SP 2 



(S41) 



In the absence of the the imaginary transverse tunneling (A y = 0), time- reversal symmetry is also respected. Here 
we have sm(6(k)) = — sin(#(— fc)) because h y (k) = —h y (—k). Apparently § dksin(9(k)) = 0, so we conclude dp 2 = | 
for A y = 0, regardless of t s , t p and t sp . 
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S-7. THE FERRO-ORBITAL ORDER OF THE MOTT INSULATOR 



The quantum phase in the strong interaction regime U sp > 11° is a Mott insulator. The double occupancy in this 
phase is found to be greatly suppressed, i.e., (n s (j)n p (j)) <C 1. The orbital physics of this Mott insulator can be 
described by an effective Hamiltonian with the double occupancy projected out. Then the two states, a\(j)\D,) and 
(— iya p (j)\il), are mapped to two pseudo-spin \ states. The resulting effective Hamiltonian (for t s = t p = t ) is the 
well-known XXZ Hamiltonian given as 



(S42) 



^2_|_^2 ^2 ^2 

with J x = 2 °jj 3p and J yz = 2 ° p sp . The XXZ Hamiltonian predicts a gapped Mott insulator for the sp-orbital 
ladder with a ferro-orbital order (C-a x Cj), which spontaneously breaks the particle-hole symmetry. 



S-8. THE STABILITY AGAINST INTER-LADDER COUPLING 



In this section the stability of the topological property of the sp-orbital ladder against the inter-ladder coupling will 
be proved. Due to the experimental setup the leading inter-ladder coupling is the coupling between the B (A) chain 
and the A (B) chain of the nearest ladder. The tight binding Hamiltonian describing such a coupled two dimensional 
system is given as 



H2D — ^ ^C^TCn+x + Cr.7 1 ' 'Cn+x+y — CrT'Cr+S-x + h.c.\ 
R 



(S43) 



where R labels the positions of lattice sites and x (y) is the primitive vector in the x- (y-) direction. The inter-ladder 

, with t' the inter-ladder coupling strength. In the momentum space the 



coupling matrix T" is given as 
Hamiltonian reads H 2 d = J2k C'^(k)'H2Li(k)C'(k), with 



t' sp 



H 2D (k) 



2ig cos kx 
2i(t sp + if e ik y) sin k x 



2i(t sp + t'e- lk v)smk 



sp 

*2itp cos 



(S44) 



The Hamiltonian "H2u(k) can be rewritten as %2.o(k) = U'(k y )'H2D(k x ,ky)U(ky) with 



T-L2D{k x -, k y ) 



cos kx ^^^sp sin k% 
c 2aAj sp sin kx ^2Ajp cos kx 



(S45) 



? sp e 



, and s(k y ) = &rg(t sp +t' sp e lk y). With k y fixed, 'H.2D{k x ,k y ) 



and U(k y ) = e i^f(fc»)/2 ) w here t sp = \t s 

defines a one dimensional system. Because the form of H2D(k x , k y ) is exactly the same as the one dimensional 
Hamiltonian describing the ,sp-orbital ladder in the limit of t' sp = 0, the Berry phase of r H.2D{k x , k y ) for any given k y 
is 7r. We thus conclude that the one dimensional topological insulator survives even with finite inter-ladder couplings. 

The existence of edge states in presence of inter-ladder coupling is verified by directly calculating the energy spectra 
of the 2D system on a cylinder geometry (an open cylinder in the cc-direction). From the energy spectra shown in 
Fig. [5] it is clear that the edge states are stable against inter -ladder couplings. 
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